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The thermodynamic Bethe ansatz (TBA) equation for an inte- 
grable spin chain related to the Lie superalgebra osp{l\2) is analyzed. 
The central charge determined by low temperature asymptotics of the 
specific heat can be expressed by the Rogers dilogarithmic function, 
and identified to be 1. Solving the TBA equation numerically, we eval- 
uate the several thermodynamic quantities. The excited state TBA 
equation is also discussed. 
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1 Introduction 

Recently, thermodynamics of quantum integrable spin chains related to Lie 
superalgebras received much attentions. In particular, several people [0, |^, 
^ |, 1^ studied thermodynamic Bethe ansatz (TBA) equations related to 
sl{r\s). On the other hand, study on osp{r\2s) case has begun only recently 
in refs. 0, Q (cf. ref. 0), in which we deal with the simplest osp{l\2) model 
TD|, |rT|. Namely, we have derived the F-system from the osp{l\2) version 



of the T-system ||T^, and transformed it into the TBA equation and excited 
state TBA equation from the point of view of the quantum transfer matrix 
(QTM) method [0. As for the largest eigenvalue sector of the dressed 
vacuum form (DVF), this TBA equation coincides with the one [0 from the 



traditional string hypothesis |T^, ^]. On the other hand, the excited state 
TBA equation from the second largest eigenvalue sector, which characterizes 
the correlation length, is difficult to derive by the string hypothesis. The 
purpose of this paper is applications of our previous results to the 

calculation of physical quantities. 

The Hamiltonian |]16| of the osp(l|2) integrable spin chain, which we deal 



with in this paper is given by 

■ ' ' (1.11 



where we assume the periodic boundary condition. For precise definition of 
( |1.1| ), see refs. [1^, ||. In this paper, we consider the case J = — 1, which 
corresponds to the antiferromagnetic regime. 

In section 2, we consider the finite temperature correction of the DVF 
T^^\v) from the T-system |12] and calculate the central charge from the 



low temperature asymptotics of the specific heat (see, refs. O, |I9|, |20, 
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2T|). The central charge is expressed in terms of the Rogers dilogarithm 
function, and reproduces the conjecture [|l6l c = 1 from the root density 
method. Moreover we solve the TBA equation numerically and evaluate the 
thermodynamic quantities such as the free energy, the internal energy, the 
specific heat and the entropy, which are depicted in fig. |l|. To the author's 
knowledge, this is the first evaluations of the physical quantities of osp{l\2) 
model at finite temperature. 

Above calculation corresponds to the largest eigenvalue sector of the DVF. 
In section 3 we comment on the excited state TBA equations, which have 
singularities from zeros of fused QTMs. For comparison, we also briefly 
mention a finite temperature correction from the traditional string hypothesis 
in Appendix. 

2 Analysis of the TBA Equation 

In our previous papers [0, ||, we proposed the TBA equation of the osp{l\2) 
model 

= + ^ * 1 TTy^^ 1 ^^-'^ 

where m G Z>i, lo(f) := 0, /? = 1/T (T: temperature; in this paper, we set 
the Boltzmann constant to 1.), K{v) is the kernel defined by 



K{v) = (2.2) 

2 cosh TTV 



and A* B(v) denotes the convolution 



oo 

A*B{v)= I dwA{v -w)B{w). (2.3) 

— oo 
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Through this TBA equation, the logarithm of the largest eigenvalue of the 
QTM Ti^\v) (See, ref.|| for the definition of t4^' (t>).) can be expressed as 

logTf )(0) = P -l]+G* log(l + ri)(0), (2.4) 



.3^3 

and the free energy per site / is written as 

/ = -Tlogrf)(0). (2.5) 

Here the kernel G{v) is defined by 

^ _ 2sinh^ ^2 
-\/3 sinh 2ttv 

Note that -47r/(3V3) + 1 =: E^^ in eq. ( p.4|) is the ground state energy 
of (O- 



Following ref. |T^, we shall analyze the TBA equation ( |2.1| ). In the 



second term of the right hand side of (|2.1|) , we find that the factor Y^^{v) is 
not relevant for the calculation of the finite temperature correction of tI^^\v). 
To avoid this, we have to modify ( p.ll ). After some manipulation, one can 
rewrite (12. 11) as 



\ogY^{v) = + G * log <^ — } (v), (2.7) 

Vosmh27rt> I, i + -f m 



where m G Z>i, lo(f) := 0. We introduce the scaling function in the low 
temperature limit 

lunK„,(±i-(„ + log M)). (2.8) 

From the modified TBA equation ( p.7|) , we find that ym,±{v) satisfies the 
following non-linear integral equation 

1 ^ N -"X , pi 1 / (1 + ?/m+l,±)(l + ym-l,±) 1 / N 

log ym,±{v) = -e S„^i + G*\og< — > [v), (2.9) 

I 1 + ym,± ) 
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where m G Z>i, yo,±{v) := and 

We shall divide tI^\v) into the ground state and the finite temperature 
correction parts, T^^{v) and rf°(f), respectively, 

Ti^\v) =Tf'{v)Tl''{v). (2.11) 

By using the functions (p.8|) , one can derive the asymptotic behavior of the 
finite temperature correction part for large f3: 

logTl'^iv) = G*log(l + Fi)(t;) 

V3 r . ^i^hf (.-^(^ + log^) 

= — / aw ) J 

^ -^-log^ sinh27rfw-|;fw + log^ 



X 



+ — / dw ) ) — 

^ J-ios^ sinh27r + fw + log^ 

= ^je^'^y d«;e-'"log(l + yi,+ H) 

du;e-'"log(l + yi,_H)| + o(-). (2.12) 

Taking account of the relation G{—v) = G{v), we can derive the following 
relation from (2.9): 

dwe"''log(l + ?/i,±(i;)) = 2 |log(l + ?/„,±(t;))— log?/m,±(^^) 
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y^-±M^^Pog(l + y) logy 



m- 

OD 



1 Jym.±{~00) 

ym,±{00 



E Ml 



L 



1 + y 

ym,±{-00) 



(2.13) 



where we assume that ym,±{v) are non-decreasing functions on f G M and 
L[v) is the Rogers dilogarithm function 

log(l-?/) logy 



L{x) 



dy 



+ 



for < X < 1. (2.14) 



Substituting ( |2.13| ) into ( |2.12| ), we obtain 



logr^'^(^) 



2w 

e 3 



m=l ^ 



m=l 



"E Ml 



L 



-oo 



1 + y.m ~{-oo, 



ym,+ {-oo) 
1 + ym,+ {-oo] 



+ o(-).(2.15) 



In our case, both ym,+ {v) and ym-{v) behave in the same manner, thus we 
set ym{v) := ym,±{v). For small T, the leading term of the specific heat C 



dT \ dT \T 



^(T=^.ogT;«,0) 



(2.16) 



is proportional |T9|, ^ to the central charge c: 

-KcT 



C 



(2.17) 



where the Fermi velocity is fp = 27r/3 1T6|. Thus, we can express the central 
charge as 



6 



m=l 



E Ml 



+ ym(oo) 



L 



1 + 2/m(-00) 



(2.18) 



Above expression is widely seen in the model related to rank one algebras 
(see for example, refs. |I7[ 0, g). 
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Now we shall evaluate the limit ym{^oo). For t> ^ oo, ( p^ ) reduces to 
the constant y-system[0 

Vm = for m G Z>i, (2.19) 

1 + ?/m 

where yo := 0. Thus, we expect (see, eq. (4.1) in ref. 0) 

,„„ = ai^IIilR. (2.20) 



The divergence from the first term in rhs of eq. (|2.9| ) in the limit v — » — oo 
is expected to be canceled by Ihs if 

yi{v) — > +0 for V -oo. (2.21) 

Then remaining ym{v) for m G Z>2 obey the constant F-system (|2.19|) . Thus 



we expect ym{—oo) is given as the solution of (|2.19| ) with m ^ m — 1: 



(fa — iVm + 2) 

lim y„,{v) = ^ ^ for m G Z>2. (2.22) 



Using these relations (|2l^ ), (IBp , (12^2^ ) and the fact L(l) = n^/G, we 



finally arrive at c = 1. This agrees with the conjecture [|T^ by the root 
density method. 



Now we shall solve the TBA equation ( p.7|) numerically to analyze the 
thermodynamic quantities at finite temperatures. To analyze this equation 
( p.7[ ) numerically, one needs to truncate it as a finite number (we call this 
number mt) of nonlinear integral equations since it is composed of an infinite 
number of nonlinear integral equations. The validity of this approximation 
can be verified from the fact that the numerical results are almost indepen- 
dent of mt when mt is more than a certain large number. 

Figure |l| is numerical results of temperature dependence of some im- 
portant thermodynamic quantities such as the free energy ( p.5|) , the spe- 
cific heat ( p.l6[ ), the entropy S = ^(TlogT^''^^(0)) and the internal energy 
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Figure 1: Temperature dependence of fundamental thermodynamic quanti- 
ties: (a) Free energy f(T), (b) Specific heat C(T), (c) Internal energy U(T) 
and (d) Entropy S(T). 



U = T^^{logTi^\0)). These results agree with well-known values at the 
both special limits: (i) at the low temperature limit, the free energy and the 



internal energy agree with the ground state energy [|T6| Eg^ ~ —1.4184 and 
the specific heat shows proportionate increase with respect to the tempera- 
ture, which leads the central charge c = 1, (ii) at the high temperature limits, 
the entropy agrees with the value log 3 ~ 1.09861 derived from the fact that 
the present model has three states. 
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3 Comment on the Excited State 

The excited state TBA equations are characterized by zeros of DVF Tm ^ (v) 
of the fused QTM and a phase factor. They have the following form (cf. ref. 

§)• 

^°s^-(")=-^;^+^*^°n — tty::^' — j^^) 



+ log 



!!{.,„} tanhf(t;-2;„ 



+{l-(-l)'"}{l-(-l)^-}-, (3.1) 

where m G Z>i; Yo{v) := 0; Coo = liniAr^oo(^ — n); N: the Trotter num- 
ber; n: a quantum number in the Bethe ansatz equation (BAE); {zm} are 
zeros of Tm\v) in the physical strip Imv G [—1/2, 1/2], which can also be 
characterized as 

YUzm±'-) = -l. (3.2) 

For the second largest eigenvalue case (n = — 1, {zm} = {^Xm\m G Z>i}: 
there is a misprint in the phase factor of eq.(5.3) in ref. |^.), we have to 
consider zeros of Tm\v), which enter into the physical strip for all m G Z>i. 
In addition, the second largest eigenvalue is conjectured to be described by 
the same distribution pattern (see. Fig. 3 and Fig.5 in ref. [§) of the root 
of the BAE at any finite temperatures. On the other hand, for the third 
largest eigenvalue case, our numerical study for finite A^ indicates that the 
distribution pattern of the roots of the BAE may vary with temperature. 
This suggests a possibility of the occurrence of the level crossing. We have 
also observed numerically that only the zeros (v) appear in the physical 

strip for one of the distribution patterns for the third largest eigenvalue. 



^This misprint was corrected in math-ph/9912014 v3 
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Numerical studies about such excited states are also a crucial problem. 
However to solve the excited state TEA equation, we have to analyze not 
only infinite number of unknown functions but also their zeros which are de- 
termined by subsidiary condition ( |3.2|) . To make numerical results converge, 
one needs to determine all the zeros in the physical strip accurately. For this 
reason the truncation approach available for solving the TEA equation have 
no use for the excited state case. 
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Appendix: Finite Temperature Correction from 
the String Hypothesis 

We can also confirm the central charge c = 1 based on the string hypothesis 



(see, refs. [^, ^ ^ ^^^^ case, the F-function in TEA ( p.7|) corre- 

sponds to the ratio of the particle density Pm(^) the hole density p^{v) 
: Ym{v) = Pm(^)/Pm(^)- From our previous results (eq. (15) in ref. 0), the 
particle and hole densities satisfy the following relation. 
26 sinh 

T. =pUv)+ptiv)+G*pliv)-G*pt__,iv)-G*pt^,iv), (A.1.1) 

V 6 Smh ZTTV 

where we define pjj = 0. To consider the low temperature behavior of ( p.7|) . 



we set Ym{v) = exp(/3em(f )) and shift the variable v — > ^{v + log For 
T -C 1, ( |2.7| ) can be expressed in terms of the function ipmiv) := /3em(^(f + 
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e'^Smi = log(l + exp(-v9„(i;))) -log(l + exp(v9„(t;))) 

-G * log(l + exp{(pm))iv) + G* log(l + exp{(pm-i))iv) 
+G*\og{l+exp{ip^+i)){v). (A.1.2) 

This is effectively the same equation as ( |2.9| ). Requiring the same shift for 
( |A.1.1|) and comparing the result with the f -derivative of ( [A.1.2|) , we find the 



asymptotics of the particle and hole densities in the limit v —* oo: 

pi{v) ^ ^(1 - f(peUv)))^eUv), (A.1.3) 



where f{<f) = (1 + e*^) ^ is the Fermi distribution function. Using ([A.1.3|) 



and taking into account the contribution from the negative large v derived 
by the shift v ^^l"^ + l^S ^'^'^ entropy per site (eq. (21) in 

ref. 0) have the following asymptotics 

, 

m=l ■^Vmi-oo) 

) 

J2 W - /(^m(oo))) - L(l - /(^„(-oo)))} . (A.1.4) 

m=l 

Using the fact C = T{dS/dT), eq. ( |2.17|) and the identification ym{v) = 



O °° /•ipm{oo) 

^ = -^E/ d^{/(y.)log/(y.) + (l-/M)log(l-/M)} 



limT^oexp(y9m('y)), we can reconfirm the central charge c = 1. 
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